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ABSTRACT 

We  have  completed  preliminary  work  to  develop  a  new  model  for  turbulent 
wall-bounded  flows  that  captures  the  correct  scaling  behavior  of  the  turbu¬ 
lence  intensities  and  other  near-wall  statistics  (including  skin  friction)  for  a 
wide  range  of  Reynolds  numbers.  Such  information  can  be  used  for  predic¬ 
tion  purposes  and  may  be  the  basis  of  a  new  approach  to  the  near-wall  model 
problem  in  Large  Eddy  Simulation  (LES). 

The  model  is  based  on  our  new  understanding  of  turbulent  structure  in 
subsonic,  supersonic  and  hypersonic  flows.  Specifically,  we  identify  four  basic 
eddy  motions:  (1)  the  streaks  associated  with  longitudinal  vortex-like  struc¬ 
tures  in  the  near-wall  region,  as  identified  by  Kline  et  al.  (1967);  (2)  the 
hairpin  or  horseshoe  vortices  first  described  by  Theodorsen  (1952);  (3)  the 
Large-Scale  Motions  (LSMs)  which  are  related  to  the  groups  of  horseshoe 
vortices  or  “vortex  packets”  defined  by  Head  &  Bandyopadhyay  (1981)  and 
Adrian  et  al.  (2000);  and  (4)  the  Very  Large-Scale  Motions  (VLSMs)  inter¬ 
preted  by  Liu  et  al.  (2001)  in  terms  of  the  outer  layer  bulges  and  by  Monty 
et  al.  (2007)  in  terms  of  the  meandering  superstructures  observed  in  pipe, 
channel  and  boundary  layers. 

We  use  the  concept  of  the  attached  eddy  hypothesis  to  map  the  attributes 
of  each  eddy  type  in  physical  space  to  wavenumber  space.  Experimental  data 
are  then  used  to  determine  the  scaling  behavior  of  the  basic  eddy  motions  in 
wavenumber  space,  and  the  scaling  behavior  of  the  Reynolds  stress  behavior 
is  recovered  by  integrating  over  all  wavenumbers. 
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1.  STATEMENT  OF  OBJECTIVES 


We  report  the  preliminary  results  of  a  combined  experimental/theoretical 
program  to  study  the  behavior  of  wall-bounded  turbulent  flows  over  a  very 
large  Reynolds  number  range.  The  aim  is  to  develop  a  new  model  for  the 
turbulence  behavior  that  is  based  on  our  knowledge  of  the  coherent  structures 
observed  in  such  flows.  The  approach  taken  here  is  to  build  on  the  attached 
eddy  model  originally  proposed  by  Townsend  (1976)  and  refined  by  Perry  & 
Chong  (1982),  Perry,  Henbest  &  Chong  (1986),  Fernando  &  Srnits  (1988), 
Perry  &  Marusic  (1995),  and  Marusic,  Uddin  &  Perry  (1997). 

The  attached  eddy  model  describes  turbulence  in  wall-bounded  flows  as 
consisting  of  a  superposition  of  a  number  of  distinct  classes  of  eddies,  where 
the  different  classes  correspond  to  the  coherent  structures  observed  in  ex¬ 
periment.  The  eddies  are  scaled  into  hierarchies,  and  each  hierarchy  has  the 
correct  scaling  and  statistical  properties  so  that  the  ensemble  of  hierarchies 
reproduces  both  the  mean  flow,  and  the  spectral  behavior  of  the  turbulence. 
The  most  important  aspect  of  the  attached  eddy  model  is  the  underlying 
scaling  argument  that  anchors  the  model  to  a  firm  physical  foundation.  Be¬ 
cause  the  model  is  based  on  eddy  physics  derived  from  observations  of  tur¬ 
bulent  structure,  it  is  expected  to  be  more  robust  to  changes  in  boundary 
conditions  than  traditional  models  used  in  Reynolds-Averaged  Navier-Stokes 
solvers.  Here,  we  build  on  this  foundation  to  define  eddy  functions  that  cap¬ 
ture  the  contributions  of  the  principal  coherent  motions  to  the  turbulence 
energy  and  shearing  stress  in  wavenumber  space.  Preliminary  results  (Srnits, 
2010)  have  demonstrated  the  remarkable  result  that  only  three  eddy  func¬ 
tions,  corresponding  to  the  large-scale  motions,  very  large-scale  motions,  and 
the  motions  responsible  for  the  near-wall  energy  production,  are  sufficient  to 
capture  the  spectral  behavior  of  wall-bounded  turbulence,  and,  by  integra¬ 
tion,  predict  the  turbulence  intensity  distribution  in  boundary  layers  over  a 
very  large  Reynolds  number  range. 

The  existing  attached  eddy  model  has  had  some  notable  successes,  such 
as  the  derivation  of  scaling  laws  for  the  three  components  of  turbulence, 
and  the  successful  scaling  of  the  streamwise  component  of  the  turbulence  for 
Reynolds  numbers  corresponding  to  small  laboratory  facilities  all  the  way 
to  atmospheric  boundary  layers  (Marusic  &  Kunkel,  2003).  It  was  extended 
to  compressible  flows,  at  least  at  an  elementary  level,  by  Fernando  &  Srnits 
(1988),  Dussauge  &  Srnits  (1995)  and  Srnits  &  Dussauge  (2006).  This  work 
has  established  the  attached  eddy  model  as  the  preferred  interpretation  for 
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the  structure  of  turbulence,  and  the  primary  basis  for  understanding  the 
interaction  among  different  scales  of  motion.  We  now  need  to  test  its  full 
capability  by  expanding  the  model  scope  so  that  it  will  reproduce  many 
other  features  of  flat  plate  boundary  layers,  such  as  the  spectral  distribution 
of  scales,  higher  order  moments,  and  components  of  turbulence  other  than 
the  streamwise  stress,  in  particular  the  shear  stress  distribution.  In  addition, 
we  now  know  that  turbulent  boundary  layers  and  other  wall-bounded  flows 
such  as  pipe  or  channel  flows  display  previously  unsuspected  differences  that 
have  not  yet  been  incorporated  in  the  attached  eddy  model.  Furthermore,  the 
model  needs  to  include  the  effects  of  pressure  gradient,  wall  curvature,  and 
heat  transfer.  Finally,  incorporating  density  variations  and  compressibility  is 
still  at  a  rudimentary  stage,  and  there  have  not  yet  been  any  serious  attempt 
to  develop  the  attached  eddy  model  to  describe  the  behavior  of  supersonic 
and  hypersonic  turbulent  boundary  layers. 

In  the  following,  we  describe  the  current  state  of  the  art  in  the  attached 
eddy  model  in  more  detail,  and  summarize  our  progress  in  expanding  the 
model  to  incorporate  more  detailed  flow  physics  based  on  a  wavelet  interpre¬ 
tation  of  eddy  structure. 

2.  RESEARCH  EFFORT 

2.1  The  Structure  of  Turbulence 

In  considering  the  structure  and  scaling  of  wall  turbulence  at  high  Reynolds 
numbers  we  start  with  the  “classical”  position  where  the  boundary  layer  is 
held  to  be  composed  of  two  principal  regions  that  follow  distinct  scalings: 
a  near- wall  region  where  viscosity  is  important,  and  the  outer  region  where 
it  is  not.  On  the  basis  of  the  mean  momentum  equation,  the  velocity  and 
length  scales  in  the  near  wall  region  are  taken  to  be  uT  =  \Jrw/p  and  v /uT, 
respectively,  where  tw  is  the  wall  stress,  p  is  the  fluid  density,  and  v  is  the  fluid 
kinematic  viscosity.  In  the  outer  region,  it  is  assumed  that  the  appropriate 
length  scale  is  the  boundary  layer  thickness  <5,  or  a  scale  related  to  5,  and 
the  velocity  scale  continues  to  be  uT,  since  uT  sets  up  the  inner  boundary 
condition  for  the  outer  flow. 

As  to  the  structure  of  turbulence,  in  the  current  view  we  identify  four 
principal  characteristic  elements  (Smits  et  al,  2011).  The  first  two,  near¬ 
wall  streaks  with  a  typical  spanwise  spacing  of  about  100 vw/uT  (Kline  et  al., 
1967),  and  hairpin  or  horseshoe  vortices  with  a  range  of  scales  starting  with  a 
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minimum  height  of  100 uw/uT  (Theodorsen,  1952),  have  been  recognized  for  a 
long  time.  More  recently,  visualizations,  numerical  studies  and  experiments 
on  wall-flows  have  revealed  the  existence  of  two  new  elements,  the  so-called 
large-  and  very  large-scale  motions. 

Large-Scale  Motions  (LSM)  are  believed  to  be  created  by  the  vortex  pack¬ 
ets  formed  when  multiple  hairpin  structures  travel  at  the  same  convective 
velocity  (Zhou  et  al.  1999,  Kim  &  Adrian  1999,  Guala  et  al.  2006,  Balaku- 
mar  &  Adrian  2007).  They  are  related  to  the  features  first  seen  by  Head  & 
Bandyopadhyay  (1981)  that  consisted  of  hairpins  with  their  heads  inclined  at 
approximately  20°  to  the  wall,  but  what  has  become  clear,  especially  through 
the  work  of  Adrian  et  al.  (2000),  is  that  these  features  are  far  more  common 
and  important  than  Head  &  Bandyopadhyay  realized  at  the  time.  A  char¬ 
acteristic  feature  of  the  LSM  is  that  the  hairpin  vortices  within  the  packet 
align  in  the  streamwise  direction  and  induce  regions  of  low  streamwise  mo¬ 
mentum  between  their  legs  (Adrian  et  al,  2000;  Brown  &  Thomas,  1977; 
Ganapathisubramani  et  al.,  2003;  Hutchins  et  al,  2005;  Tomkins  &  Adrian, 
2003).  The  LSM  have  a  streamwise  scale  of  approximately  2-3<5  and  have 
been  associated  with  the  occurrence  of  bulges  of  turbulent  fluid  at  the  edge 
of  the  wall  layer,  where  6  is  the  boundary  layer  thickness.  A  thorough  re¬ 
view  of  the  evidence  supporting  the  existence  of  hairpin  vortices  and  their 
organization  into  packets  is  provided  by  Adrian  (2007). 

In  addition  to  the  LSM,  very  long,  meandering,  features  consisting  of 
narrow  regions  of  low  streamwise  momentum  fluid  flanked  by  regions  of  higher 
momentum  fluid  have  been  observed  in  the  logarithmic  and  wake  regions 
of  wall-flows  (Balakumar  &  Adrian,  2007;  Guala  et  al,  2006;  Hutchins  & 
Marusic,  2007 6;  Kim  &  Adrian,  1999;  Monty  et  al,  2007;  Tomkins  &  Adrian, 
2005).  In  internal  flows,  the  motions  are  typically  referred  to  as  Very- Large- 
Scale  Motions  (VLSM),  whereas  in  external  flows  they  are  more  commonly 
referred  to  as  “superstructures.”  Both  VLSM  and  superstructures  appear  to 
scale  on  outer  variables,  and  although  the  spanwise/ azimuthal  meandering  of 
these  regions  makes  it  difficult  to  determine  their  typical  streamwise  extent, 
hot-wire  rake  measurements  in  channels  and  pipes  (Monty  et  al,  2007)  have 
found  instances  of  VLSM  in  internal  flows  as  long  as  30  times  the  channel 
half-height  h  or  pipe  radius  R,  while  similar  experiments  in  boundary  layers 
(Hutchins  &  Marusic,  20076)  show  instances  of  superstructures  with  lengths 
up  to  10-15  times  the  boundary  layer  thickness.  These  lengths  are  typically 
shorter  when  inferred  from  single-point  frequency  spectra  (~10-20i?  for  pipes 
and  ~  6 5  for  boundary  layers).  In  addition,  Monty  et  al.  (2009)  note  that  the 
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superstructures  in  boundary  layers  appear  to  be  limited  to  the  logarithmic 
region,  whereas  for  internal  geometries  the  VLSM  are  found  to  persist  well 
into  the  outer  layer  (Bailey  &  Smits,  2010).  In  contrast,  Tntknn  et  al.  (2009) 
found  evidence  of  weak  elongated  structures  within  boundary  layers  out  to 
the  edge  of  the  layer. 

Spectral  analysis  of  the  VLSM  and  LSM  indicates  that  they  make  a  sig¬ 
nificant  contribution  to  the  turbulent  kinetic  energy  and  Reynolds  stress  pro¬ 
duction  (Balakumar  &  Adrian,  2007;  Guala  et  al.,  2006),  which  distinguishes 
them  from  the  inactive  motions  proposed  by  Townsend  (1976).  For  exam¬ 
ple,  Balakumar  &  Adrian  (2007)  found  that  40-65%  of  the  kinetic  energy 
and  30-50%  of  the  Reynolds  shear  stress  is  accounted  for  in  the  long  modes 
with  streamwise  wavelengths  A x/5  >  3.  Similar  estimates  for  the  contribu¬ 
tion  to  the  Reynolds  shear  stress  from  vortex  packet  structures  were  made  by 
Ganapathisubramani  et  al.  (2003),  who  used  a  feature-detection  algorithm  on 
stereoscopic  PIV  data  in  streamwise-spanwise  planes  in  the  log  layer.  Careful 
analysis  of  DNS  data  has  also  revealed  the  footprint  of  the  outer-scaled  mo¬ 
tions  within  the  inner-scaled  inner  layer  (Abe  et  al. ,  2004;  Hoyas  &  Jimenez, 
2006;  Hutchins  &  Marusic,  2007a, b),  and  Mathis  et  al.  (2009)  found  modula¬ 
tion  of  the  near  wall  cycle  by  long  wavelength  motions  further  from  the  wall, 
which  is  supported  by  the  correlations  measured  by  Tntknn  et  al.  (2009)  and 
the  recent  simulation  results  by  Schlatter  et  al.  (2009).  These  results  suggest 
that  motions  in  the  logarithmic  and  outer  layer  may  have  a  strong  influence 
on  the  behavior  of  the  near  wall  turbulence.  This  is  a  significant  observation 
in  that  it  suggests  a  new  approach  for  developing  near- wall  models  for  LES, 
a  topic  that  we  will  explore  as  an  important  part  of  our  future  work. 

2.2  The  Attached  Eddy  Model 

Here  we  give  a  description  of  the  attached  eddy  model,  and  how  the  obser¬ 
vations  on  coherent  motions  are  incorporated  in  the  model.  The  model  es¬ 
sentially  provides  a  kinematic  description  for  wall-bounded  turbulence.  The 
foundation  of  the  model,  which  is  based  on  the  attached-eddy  hypothesis  of 
Townsend  (1976),  is  the  observation  that  wall-bounded  turbulence  contains 
a  collection  of  coherent  structures  or  eddies.  Therefore,  the  model  proposes 
that  the  statistical  features  of  wall-bounded  turbulence  can  be  modelled  by 
a  linear  superposition  of  such  eddies.  The  model  has  been  refined  and  devel¬ 
oped  over  the  past  three  decades  on  the  data  of  low-to-moderate  Reynolds 
number  wall-bounded  turbulence  experiments  (e.g.  pipe  and  boundary- layer 
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flow)  and  has  led  to  a  number  of  similarity  laws.  For  a  complete  review  of  the 
attached-eddy  model  see  Perry  &  Marusic  (1995)  and  the  references  therein. 

One  of  the  strengths  of  the  attached-eddy  model  is  the  physical  framework 
it  provides  for  the  analysis  and  understanding  of  the  mechanism  of  wall- 
bounded  turbulence.  For  example,  Perry  &  Abell  (1975),  Perry,  Henbest  & 
Chong  (1986),  and  Perry  &  Li  (1990)  used  this  framework  along  with  simple 
dimensional-analysis  arguments  applied  to  the  velocity  contributions  of  the 
attached  eddies  to  suggest  several  regions  of  spectral  scaling.  They  assumed 
that  a  turbulent  boundary  layer  outside  the  viscous  region  may  be  modeled 
as  a  forest  of  hairpin  or  A-shaped  vortices,  which  originate  at  the  wall  and 
grow  outward.  Figure  1  shows  three  A-shaped  vortices  of  different  scales, 
and  indicates  their  influence  on  the  velocity  field  sensed  by  a  probe  at  a 
position  y.  The  probe  will  sense  contributions  to  u'  and  w'  from  all  eddies 
of  scale  y  and  larger.  However,  only  eddies  of  scale  y  will  contribute  to  v' . 
Therefore,  u'  and  w'  should  follow  similar  scaling  laws,  whereas  v'  is  expected 
to  follow  a  somewhat  different  scaling  law.  LIsing  these  ideas  in  conjunction 
with  dimensional  analysis,  scaling  laws  can  be  derived  for  the  energy  spectra 
in  the  turbulent  wall  region,  defined  as  v /uT  -C  y  <C  5  (Perry  et  ah,  1986).  In 
general,  it  is  the  region  where  it  is  assumed  that  direct  wall  effects  such  as  the 
damping  of  the  velocity  components  are  unimportant,  and  where  the  direct 
influence  of  the  large-scale  flow  geometry  and  outer  boundary  conditions  can 
also  be  neglected.  This  region  corresponds  approximately  to  the  overlap 
region  in  the  mean  velocity  profile,  that  is,  the  region  where  the  velocity  has 
a  logarithmic  variation. 

For  the  '(//-component  of  the  turbulence  fluctuations  in  an  incompressible 
flow,  it  is  argued  that  eddies  of  scale  5  will  contribute  only  to  the  large-scale, 
low  wavenumber  (low  frequency)  region  of  the  energy  spectrum  $n.  For 
the  large-scale  eddies,  viscosity  is  not  important,  and  the  spectrum  in  the 
low  wavenumber  region  should  depend  only  on  uT,  k±,  y ,  and  <5,  where  k\  is 
the  streamwise  component  of  the  three-dimensional  wavenumber  vector  k. 
Eddies  of  scale  y  (the  attached  eddies)  will  contribute  to  the  intermediate 
wavenumber  range  of  the  spectrum,  whereas  eddies  of  scale  5  will  not.  The 
smallest  scale  motions,  which  contribute  to  the  high  wavenumber  range  of  the 
spectrum,  depend  on  viscosity.  Kolmogorov  (1941)  assumed  that  these  small- 
scale  motions  are  locally  isotropic,  and  that  their  energy  content  depends 
only  on  the  local  rate  of  turbulence  energy  dissipation  £,  and  the  kinematic 
viscosity  z/,  so  that  their  characteristic  length  scale  is  77  =  (z/3/^)1/4. 

Just  as  the  mean  flow  exhibits  an  inner  and  outer  scaling  with  a  region 
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Figure  1:  Sketch  of  the  streamline  patterns  and  spatial  influence  of  attached  eddies 
at  three  different  scales.  From  Perry  et  al.  (1986). 


of  overlap,  it  is  expected  that  there  will  be  an  overlap  region  in  wavenumber 
space  between  the  5-sized  motions  and  the  //-sized  motions.  Dimensional 
analysis  then  indicates  that,  in  this  overlap  region,  the  spectrum  follows  a  k~l 
scaling.  We  also  expect  that  there  will  be  an  overlap  region  in  wavenumber 
space  between  the  //-sized  motions  and  the  //-sized  motions.  This  leads  to  a 
spectrum  following  a  /c-5/3  scaling. 

By  integrating  these  spectral  forms,  Perry  &  Li  (1990)  derived  the  fol¬ 
lowing  expressions  for  the  normal  stresses  in  a  subsonic  boundary  layer: 

%  =  -  A,  !n  (|)  -  V(y+),  (1) 

|  -  *-*-(3)-™  (2) 

V4  =  A2~V(y+),  (3) 

where  B\  and  B2  are  large-scale  characteristic  constants,  particular  to  the 
flow  geometry,  and  A\,  A2,  and  A3  are  expected  to  be  universal  constants  (all 
the  constants  are  positive).  The  logarithmic  term  comes  from  the  A;-1  portion 
of  the  spectrum.  The  function  V(y+ )  comes  from  the  inertial  (. k~ 5//3  plus 
dissipation)  portion  of  the  spectrum,  and  it  is  a  Reynolds  number  dependent, 
viscous  correction  term  that  increases  with  y+ .  Note  that  equations  1  to  3 
are  valid  only  in  the  turbulent  wall  region. 
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These  turbulence  intensity  similarity  formulations,  and  the  spectral  sim¬ 
ilarity  arguments  on  which  they  are  based,  have  been  tested  for  a  variety  of 
wall-bounded  flows  and  the  results  used  to  further  improve  the  model.  For  in¬ 
stance,  the  suggested  formulations  have  been  investigated  using  smooth-pipe 
flow  (Perry  et  al.  1986),  smooth  and  rough  zero-pressure-gradient  turbulent 
boundary  layers  (Perry  &  Li  1990),  and  favorable-  and  adverse-pressure- 
gradient  turbulent  boundary  layers  (Jones,  Marusic  &  Perry  2001  and  Maru- 
sic  &  Perry  1995,  respectively). 

Marusic  et  al.  (1997)  extended  these  scaling  laws  to  include  the  entire 
region  outside  the  viscous  sublayer  (the  original  formulations  were  only  valid 
in  the  log  region).  As  with  the  mean  flow,  the  deviation  from  the  logarithmic 
profile  near  the  wall  is  attributed  to  viscous  effects,  and  the  deviation  in  the 
outer  part  of  the  layer  is  due  to  wake  effects.  They  suggested  a  “wall- wake” 
distribution  where,  for  example, 

(4) 

They  called  Vg\  the  viscous  deviation  (actually  a  more  complete  version  of 
function  V  in  equation  1),  and  Wgl  the  wake  deviation.  Marusic  et  al.  (1997) 
gave  empirical  forms  for  Vgi  and  Wg i  that  agreed  well  with  data  over  the 
range  6570  <  Ree  <  35100.  An  extension  of  Equation  4  that  applies  across 
the  entire  boundary  layer  including  the  viscous  near-wall  layer  was  proposed 
by  Marusic  &  Kunkel  (2003),  and  comparisons  with  the  laboratory  data  of 
DeGraaff  &  Eaton  (2000)  and  the  atmospheric  data  of  Metzger  &  Klewicki 
(2001)  showed  excellent  agreement  (Figure  2),  although  the  formulation  does 
not  reproduce  the  outer  peak  in  the  streamwise  intensity  observed  at  high 
Reynolds  numbers  seen  by  Fernholz  et  al.  (1995),  and  Morrison  et  al.  (2004). 
A  similar  high  level  of  agreement  with  high  Reynolds  number  pipe  flow  data 
was  shown  by  Marusic  et  al.  (2004). 

Despite  its  apparent  success,  there  are  a  number  of  compelling  reasons 
to  revisit  this  model.  First,  the  model  is  based  on  the  presence  of  /c_1  and 
k~5/3  regions  in  the  spectrum.  Although  the  fc-5/3  is  well  established,  experi¬ 
ments  now  indicate  that  the  /c-1  region  is  only  evident  at  very  high  Reynolds 
numbers  over  a  very  limited  spatial  extent  (see,  for  example,  Morrison  et 
al.  2004,  and  Nickels  et  al.  2005).  Second,  the  interactions  between  outer 
layer  motions  and  inner  layer  motions  have  become  much  clearer  in  recent 
years,  and  the  simple  division  between  inner  and  outer  layer  scaling  that 


Figure  2:  Streamwise  turbulence  intensity  measurements:  high  Reynolds  numbers 
from  ASL  (Metzger  &  Klewicki;  Marusic  et  al.);  low  Reynolds  number  data  from 
DeGraaff  &  Eaton  (2000).  Lines  are  model  results  from  Marusic  &  Kunkel  (2003). 

leads  to  the  k  l  region  fails  to  capture  that  interaction.  Specifically,  the 
region  where  we  might  expect  kl  scaling  corresponds  to  the  wavenumbers 
occupied  by  the  LSMs,  and  experiments  have  clearly  shown  that  although 
the  LSMs  appear  to  behave  as  attached  motions,  they  do  not  scale  simply  as 
y -1.  Third,  the  importance  of  the  VLSMs  was  not  appreciated  until  recently. 
For  example,  the  low  wavenumber  VLSMs  contribute  about  half  of  the  total 
energy  content  of  the  streamwise  turbulence  component  at  high  Reynolds 
numbers.  Fourth,  it  has  become  clear  that  the  relative  importance  of  LSMs, 
VLSMs,  and  superstructures  depends  on  the  nature  of  the  flow:  they  behave 
differently  in  pipes,  channels  and  boundary  layers  (Monty  et  al.  2007,  Bailey 
et  al.  2008),  something  that  is  not  captured  in  the  original  model.  Fifth, 
the  original  model  does  not  include  the  effects  of  pressure  gradient  (although 
some  steps  were  taken  by  Perry,  Marusic  &  Jones  2002),  compressibility,  or 
heat  transfer. 

As  to  the  effects  of  compressibility,  recent  results  at  Princeton  on  the 
behavior  of  hypersonic  boundary  layers  have  given  new  confidence  that  such 
flows  can  be  incorporated  in  the  attached  eddy  approach.  We  have  demon- 
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Figure  3:  Intensity  of  the  streamwise  component  of  the  velocity  fluctuation  in  a 
Mach  7.2  turbulent  boundary  layer  using  Morkovin  scaling  (Rer  =  300,  Reg  = 
3500).  PIV  data  from  Sahoo  et  al.  (2010a),  DNS  from  Duan  et  al.  (2010). 

strated  that  the  mean  velocity  profile  for  smooth  and  rough  surfaces,  when 
scaled  according  to  van  Driest,  follows  the  incompressible  scaling  closely  (Sa¬ 
hoo  et  al.,  20106,  2009).  In  addition,  the  experimental  data  of  Sahoo  et  al. 
(2010a)  and  the  DNS  results  by  Duan  et  al.  (2010)  have  provided  the  first 
clear  evidence  that  Morkovin’s  scaling  continues  to  apply  up  to  Mach  num¬ 
bers  of  7.2  (see  Figure  3).  These  results  validate  the  approach  taken  by 
Fernando  &  Smits  (1988)  in  extending  the  attached  eddy  hypothesis  to  com¬ 
pressible  boundary  layers,  and  Dussauge  &  Smits  (1995)  in  extending  the 
scaling  arguments  underlying  equation  1  to  high  speed  flows.  Taken  as  a 
whole,  this  work  opens  up  a  clear  path  to  generalizing  the  incompressible 
model  to  high  Mach  numbers. 

2.3  Technical  Approach  and  Results 

To  overcome  the  limitations  of  the  existing  model  described  above,  we  have 
conducted  a  combined  theoretical  and  experimental  study  to  develop  a  new 
model  that  captures  the  correct  spectral  behavior  of  wall-bounded  turbu¬ 
lence,  and  allows  prediction  of  Reynolds  stress  distributions  in  general  wall- 
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Figure  4:  Contour  maps  of  pre- multiplied  spectra  for  the  streamwise  component 
of  the  velocity  fluctuation  in  a  zero  pressure  gradient  turbulent  boundary  layer  for 
ReT  =  14K.  Left:  data  from  Hutchins  &  Marusic  (20076);  right:  model  predictions. 

bounded  flows. 

A  preliminary  form  for  the  new  model  has  now  been  established  (Smits, 
2010).  Here,  the  four  basic  eddy  motions  are  modeled  using  Gaussian  dis¬ 
tributions  in  the  log  of  the  wave  number.  This  represents  a  more  complete 
formulation,  when  compared  to  the  Perry  et  al.  model,  in  the  sense  that  the 
full  wavenumber  distributions  are  specified,  instead  of  just  noting  regions  of 
k~l  and  k~ 5'3  scaling.  The  energy  contained  at  each  wavenumber  is  taken 
to  be  a  simple  sum  over  all  eddy  types,  neglecting  nonlinear  interactions 
among  eddies.  The  high  wavenumber  part  of  the  spectrum  is  bounded  by 
the  fc-5//3  inertial  and  dissipation  ranges,  which  are  modeled  as  a  single  cutoff 
using  a  modified  Pao  spectrum.  The  results  are  compared  with  the  data  of 
(Hutchins  &  Marusic,  20076)  in  figure  4,  which  show  the  contour  maps  of 
pre-multiplied  spectra  for  the  streamwise  component  of  the  velocity  fluctu¬ 
ation  in  a  zero  pressure  gradient  turbulent  boundary  layer.  The  predictions 
for  the  turbulence  intensity  for  three  different  Reynolds  numbers  are  given 
in  figure  5,  which  shows  that  even  at  this  preliminary  stage  the  new  model 
performs  as  well  as  the  original  model  shown  in  figure  2  while  capturing  the 
actual  spectral  behavior  much  more  closely. 

The  results  are  remarkable  in  that  only  three  eddy  functions  are  sufficient 
to  capture  the  spectral  behavior  of  wall-bounded  turbulence,  and,  by  inte¬ 
gration,  predict  the  turbulence  intensity  distribution  in  boundary  layers  over 
a  very  large  Reynolds  number  range.  The  model  shown  in  figure  4  uses  sim- 
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Figure  5:  Intensity  of  the  streamwise  component  of  the  velocity  fluctuation.  Data 
points  from  Bailey  (pipe,  Rer  =  3.3K),  DeGraaff  (boundary  layer,  Rer  =  14K), 
and  Kunkel  (boundary  layer,  ReT  =  2.3M),  respectively.  Lines  are  model  predic¬ 
tions  at  the  same  Reynolds  numbers. 


pie  Gaussian  distributions  in  log-wavenumber  space  to  represent  the  energy 
content  of  the  three  basic  eddy  structures  (near-wall  streaks,  hairpin  vortices 
arranged  in  packets  to  form  LSMs,  and  VLSMs/superstructures),  and  em¬ 
pirically  determined  functions  are  used  to  control  the  scaling  of  these  energy 
distributions  in  physical  and  wavenumber  space.  The  log-wavenumber  dis¬ 
tributions  used  here  have  some  attractive  properties  (they  are  analytically 
integrable,  for  example),  but  we  need  to  develop  eddy  functions  that  have  a 
more  physical  basis.  We  are  attracted  to  the  concept  of  wavelet-based  eddy 
functions,  primarily  because  wavelets  appear  to  represent  well  the  character¬ 
istic  velocity  signatures  seen  in  turbulence.  For  example,  Tennekes  &  Lurnley 
(1972)  proposed  the  concept  of  a  “simple  eddy”  as  a  possible  model  for  the 
coherent  motions  in  turbulent  flows,  and  the  simple  eddy  they  suggested  is 
well-represented  by  a  Mexican  Hat  wavelet  function  in  physical  space,  as  we 
shall  see. 

Some  guidance  on  the  use  of  wavelets  to  build  eddy  functions  can  be 
found  in  recent  experimental  results  from  the  Princeton  University  Superpipe 
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(Bailey  &  Smits,  2010).  These  experiments  used  filtered  cross-correlations  of 
multi-probe  hot-wire  data  to  elucidate  more  information  about  the  structure 
of  LSMs  and  VLSMs.  These  filtered  cross-correlations,  samples  of  which 
are  shown  in  figure  7,  suggest  that  these  motions  could  be  modeled  using 
simple  wavelet  eddy  models.  Figure  7  indicates  that  the  form  of  wavelet  best 
representing  each  type  of  motion  varies,  with  VLSMs  represented  using 

g(t)  =  -e-d/^“)2  (5) 

a 

and  LSMs  represented  using  a  Mexican  Hat  wavelet  such  as 

l--Ve-^a)2.  (6) 

a  ) 


where  a  is  a  scaling  value. 

These  functions  can  be  considered  the  real  part  of  a  complex  function, 
g(t),  whose  Fourier  transforms  are  given  by 


G(u)  =  j(au)e-iau)2/2 


(7) 


and 


G{u) 


2(aca)2e-(aaj)2/2  for  cu  >  0 
0  for  u  <  0 


(8) 


The  wavelet  functions  g(t)  are  shown  in  figure  6  multiplied  by  a  weighting 
function  w(a)  used  to  ensure  that  the  energy  of  the  wavelet  is  the  same  at 
all  scales.  Also  shown  in  figure  6  are  the  modulus  of  the  corresponding 
Fourier  transforms.  These  representations  illustrate  the  localized  nature  of 
the  wavelet  function  in  both  physical  and  frequency  space. 

The  Fourier  transform  of  these  functions  were  used  to  model  the  eddy 
functions  instead  of  the  original  Gaussian  distributions.  As  shown  in  figure  8, 
the  eddy  distribution  representing  the  VLSM/superstructure  contribution 
was  very  successfully  modeled  by  the  wavelet  function  of  equation  5,  and 
actually  eliminated  the  need  for  one  of  the  empirical  functions  used  to  define 
the  original  Gaussian  distribution  of  energy. 

However,  our  attempts  at  modeling  the  LSM  contribution  using  a  single 
wavelet  function  were  less  successful.  As  shown  in  figure  9a-b,  a  single  wavelet 
does  not  incorporate  energy  across  a  large  enough  wavenumber  range  to 
model  the  energy  content  of  the  LSMs  in  a  simple  way.  Note  in  particular 
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(b) 


Figure  6:  Representative  wavelet  function  for  three  different  scale  factors  a  and 
w(a)  =  1/y/E:  (a)  equation  5  in  physical  space  (b)  equation  5  in  frequency  space 
(c)  Mexican  Hat  in  physical  space  and  (d)  Mexican  Hat  in  frequency  space. 

that  the  Fourier  transform  of  the  Mexican  Hat  wavelet  function  is  Gaussian 
in  the  wavenumber,  not  the  log  of  the  wavenumber,  so  that  the  eddy  functions 
used  in  future  attempts  to  model  the  spectral  behavior  of  the  eddies  may  need 
to  be  modified  accordingly.  Alternatively,  this  obstacle  can  be  overcome  if  we 
model  the  LSMs  as  a  hierarchy  of  eddies  of  containing  different  wavenumbers, 
as  is  done  with  the  Perry  et  al.  model.  A  simple  demonstration  is  provided  in 
figure  9c  where  the  probability  density  function  of  the  eddy  scales  is  modeled 
using  a  Gaussian  distribution.  These  results  illustrate  that  we  can  increase 
the  wavenumber  distribution  of  the  model  LSM  eddies,  but  this  is  done 
through  the  introduction  of  a  new  parameter  which  must  also  be  modeled: 
the  probability  distribution  of  eddies,  which  corresponds  (at  least  in  spirit) 
to  the  concept  of  a  hierarchy  of  eddies  as  used  in  the  refined  attached  eddy 
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Figure  7:  Spectrally  filtered  cross-correlations  with  reference  probe  at  different 
locations  from  the  wall  for  (a)  VLSM  wavenumber  range  and  (b)  LSM  wavenumber 
range. 


hypothesis. 

The  modeling  of  the  VSLMs  (or  superstructures),  together  with  the  re¬ 
cent  finding  of  Hutchins  &  Marusic  (2007)  documenting  the  interaction  of 
the  large-eddy  motions  on  the  near  wall  region,  have  led  to  us  to  propose 
a  new  preliminary  predictive  model  for  boundary  layers.  Here,  a  statistical 
signature  of  the  full  near-wall  region  can  be  predicted  given  only  the  low-pass 
filtered  velocity  signature  in  the  log-region  (such  as  that  found  in  LES).  The 
model  is  non-linear  and  is  based  on  a  Townsend  attached  signature  together 
with  an  amplitude  modulation  effect.  Preliminary  tests  show  the  predic¬ 
tive  model  to  be  very  promising.  For  example,  in  figure  10a  premultiplied 
streamwise  velocity  spectra  at  y+  =  15  (nominally  the  location  of  maximum 
turbulent  energy  production)  is  shown  in  a  high  Reynolds  number  experi¬ 
ment  together  with  a  prediction  using  the  new  model  based  on  the  velocity 
signature  measured  at  another  time  in  the  logarithmic  region  of  the  bound¬ 
ary  layer.  Figure  10b  shows  the  predicted  turbulence  intensities  at  y+  =  15 
for  a  range  of  Reynolds  numbers.  The  model  also  accurately  predicts  the 
increase  in  skewness  and  flatness  of  the  streamwise  velocity  fluctuations  with 
increasing  Reynolds  number  at  y+  =  15  (a  fact  that  has  been  known  for  some 
time  but  never  understood). 

Although  the  current  work  has  indicated  the  potential  of  this  model,  there 
are  still  several  areas  which  require  further  exploration: 

•  What  probability  distributions  can  be  best  used  to  describe  the  arrange¬ 
ment  of  the  eddies?  These  distributions  will  be  a  function  of  both  streamwise 


15 


(b) 


log  (y/R) 


0.6 

0.5 


0.4 


Figure  8:  Model  energy  distribution  for  VLSMs/superstructures  using  (a)  Gaus¬ 
sian  in  log-wavenunrber  space  (b)  wavelet  eddy  representation. 


and  wall-normal  wavenumbers  and  need  to  be  determined  from  experiment. 

•  Can  we  find  an  analytical  description  with  which  to  replace  the  empirical 
functions  used  to  describe  the  scaling  of  the  eddy  functions? 

•  Can  the  differences  between  the  VLSMs  and  superstructures,  and  by  ex¬ 
tension  internal  and  external  flows,  be  accurately  captured  by  a  single  model 
or  will  there  be  a  need  for  separate  models  for  these  different  flows? 

•  It  has  been  suggested  that  the  LSMs  and  VLSMs  are  inter-connected.  If 
so,  it  is  possible  that  the  simple  linear  superposition  of  energy  used  in  this 
model  could  be  a  poor  representation  of  the  physics  of  their  motions.  What 
is  the  nature  of  this  interaction  and  can  we  model  it?  If  this  interconnectivity 
can  be  determined  and  incorporated  into  the  model,  it  could  greatly  reduce 
the  number  of  parameters  used  to  describe  the  eddy  functions.  Support  for 
this  inter-connectivity  of  the  LSMs  and  VLSMs,  along  with  some  suggestions 
for  the  nature  of  their  interaction,  has  been  found  in  the  experiments  of  Bai¬ 
ley  &  Smits  (2010). 

•  The  energy  within  the  near  wall  streaks  has  been  found  to  be  modulated 
by  the  superstructures  (Hutchins  &  Marusic,  20076).  Can  this  interaction 
be  incorporated  into  the  eddy  function  describing  the  contribution  of  the 
near  wall  streaks?  Can  we  develop  this  approach  into  a  new  wall  function 
formulation  for  LES? 

•  Can  the  model  be  expanded  to  describe  the  additional  Reynolds  stress 
components,  especially  the  shear  stress? 

•  Can  the  model  be  expanded  to  incorporate  various  effects;  that  is,  wall- 
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Figure  9:  Model  energy  distribution  for  LSMs  using  (a)  Gaussian  in  log- 
wavenumber  space  (b)  single  wavelet  eddy  representation  (c)  Gaussian  distribution 
of  wavelet  model  eddies. 

roughness,  pressure  gradients,  thermal  stability? 


2.4  Synergistic  Activities 

The  work  was  performed  jointly  by  investigators  at  Princeton  University 
and  the  University  of  Melbourne,  Australia.  The  PI  has  a  long  history  of 
collaboration  with  Ivan  Marusic,  the  Melbourne  University  group  headed 
by  Professor  Ivan  Marusic,  as  co-organizers  of  the  International  Workshops 
on  Wall-Bounded  Turbulent  Flows  (2008,  2006,  2005,  2004,  and  2003),  co¬ 
participants  in  ICET  (International  Collaboration  on  Experiments  in  Turbu¬ 
lence),  and  most  recently  as  co-authors  (with  Beverley  McKeon  of  Caltech) 


17 


(IL 


8  8 


ReT 


Figure  10:  (a)  Premultiplied  spectra  of  the  streamwise  component  of  the  velocity 
fluctuation  in  a  zero  pressure  gradient  turbulent  boundary  layer  at  high  Reynolds 
number,  together  with  prediction  from  new  model  that  only  uses  the  velocity 
signature  in  the  log  region,  (b)  Model  predictions  of  turbulence  intensities  at 
y+  =  15  compared  with  data  for  range  of  Rer  (up  to  values  obtained  on  the  Utah 
salt- flats). 

of  an  article  for  Annual  Reviews  on  high  Reynolds  number  wall-bounded  tur¬ 
bulent  flows  (Smits  et  al.  2011).  Furthermore,  the  PI  has  an  appointment 
as  Honorary  Fellow  at  the  University  of  Melbourne. 
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